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1.0 INTRODUCTION 

TRW has been actively involved during the last few years in the develop- 
ment of variable conductance heat pipes for application in advanced space 
thermal control systems. In 1975, under Contract NAS2-8310 to Ames Research 
Center, TRW fabricated and tested two prototype vapor-modulated heat pipes 
which achieved variable conductance by induced wick/groove dry-out. The first 
prototype was designed for moderate capacity in a double-heat-pipe configura- 
tion. Tests of this prototype showed the roundness of the new induced-dry-out 
mechanism and uncovered some secondary elements that adversely affected per- 
formance. Recommendations for an improved design were subsequently incorporated 
in the fabrication of a second high-capacity vapor-modulated prototype. The 
new design, which uses a short vapor-modulated heat pipe to couple two conven- 
tional heat pipes, was tested successfully. The heat pipe operated at twice 
the 100-watt design heat load, and the source temperature was practically inde- 
pendent of sink conditions and increased with heat load at a rate of only 0.03°C 
per watt. The results of the preliminary tests, as well as details of the con- 
figuration of the heat pipe, are presented in a research report, CR-137782, 
prepared for Ames Research Center. 

In 1977, under Contract NAS 2-9834, TRW undertook a program to use the 
aforementioned high-capacity vapor-modulated heat pipe to study protection from 
freezing-point failure, increased control sensitivity, and transient behavior 
under a wide range of operating conditions in order to determine the full per- 
formance potential of the vapor-modulated heat pipe. In addition, a new con- 
cept, based on the vapor-induced-dry-out principle, was developed for passive 
feed back temperature control as a heat pipe diode. 

This report documents the work performed under the contract and describes 
a) the experimental and theoretical investigation of the performance of the 
vapor-modulated heat pipe, and b) the design, fabrication and test of the heat 
pipe diode. 
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2.0 


FURTHER RESEARCH ON THE VAPOR MODULATED HEAT PIPE (VMHP), 


2.1 Design. 

The details of the VMHP are presented in Research Report CR-1 37782, pre- 
pared for Ames Research Center. However, for completeness, certain basic 
design features of the heat pipe are briefly described in what follows. The 
high-capacity design has a triple-heat-pipe configuration. The tested config- 
uration is shown in Figure 2.1(c). Two conventional 1.27-cm OD, uniform poros- 
ity metal fiber slab wick heat pipes are coupled by a 2.54-cm OD vapor modu- 
lated heat pipe 15,2-cm long. The coupler heat pipe uses a wick assembly with 
a two-step gradation in porosity. The wick porosity on the condenser side of 
the bulkhead, which supports the bellows/valve assembly, is 86%, and on the 
evaporator side is 76%. Ammonia is the working fluid in all three heat pipes. 

The bellows and sensor volumes are filled with perflamopentane. 

2.2 Experimental Investigation of Transient Response of VMHP. 

2.2.1 Test Procedure, 

The heat pipe was instrumented with eighteen copper-constantan thermo- 
couples located along the pipe as shown in $K 75061 in Appendix A. All the 
temperatures, except TC 7, were monitored with a 24-channel strip chart recorder. 
Owing to the rather slow response of this instrument (1.5 to 2 minutes print- 
ing interval), a fast response single-pen millivolt strip chart recorder was 
used to continuously monitor the transient behavior of the input heat pipe 
adiabatic temperature, TC 7, Certain features of this instrument, such as its 
adjustable zero and variable span, were necessary in order to observe with high 
resolution the details of the temperature response profiles that would enable a 
better understanding of the characteristics of the vapor-induced dry out mecha- 
nism, A sample temperature profile recorded with this instrument during one of 
the test runs is shown in Figure 2.2. 

Two heat source configurations were used to simulate variations in the 
thermal mass of the evaporator section. A 1.3 Kg aluminum block with tape 
heaters attached to it, and the sensor inserted into a hole drilled into it, 
was used for the high thermal mass simulation. To simulate a thermal resistance 
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Figure 2.1. Vapor Modulated Heat Pipe Based on Induced Wick/Groove Dry Out 
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between a heat pipe and the equipment whose temperature is being controlled, 
.006 inch of Teflon was placed between the heater block and the evaporator 
saddle. 


For the low thermal mass simulation, the aluminum block was removed, and 
the tape heaters were .placed directly on top of the evaporator saddle. The 
sensor was inserted into a copper flange and mounted on the saddle between the 
tape heaters. The copper flange was designed to provide a nearly uniform tem- 
perature distribution around the perimeter of the sensor tube. The thermal 
resistance between the source and the heat pipe resulted from .006 inch of 
Teflon placed between the heat pipe and the evaporator saddle. 

The condenser section of the output heat pipe was clamped between two 
0.5-inch thick aluminum saddles and mounted on a temperature-controlled cold 
plate. To simulate a condenser thermally close-coupled to a sink, a thin film 
of RTV was placed between the condenser saddle and the cold plate. The opera- 
tion of the heat pipe in space was simulated by inserting 3/16 inch of Teflon 
between the condenser saddle and the cold plate. 

Q 2.2.2 Test Parameters. 

Several operating parameters such as the source thermal mass, condenser 
thermal resistance, sink temperature, sensor volume and power input were varied 
to determine their effects on the transient behavior of the vapor modulated 
heat pipe. The range of these parameters is shown in Table 2.1 


Table 2.1. Range of Test Parameters 


Evaporator Block Mass (M) 

.13 - 

1.3 Kg 

Control Sensor Volume (V) 

2.45 - 

4.90 cc 

Input Power 

3 - 

100 watts 

Condenser Section Thermal 
Resistance (R) 

.10 - 

.63 °C/watt 

Sink Temperature 

-57 - 

-18 °C 
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2.2.3 Test Results. 


The results from forty-two test runs are summarized in Table 2.2. Selected 
transient profiles of the input heat pipe adiabatic temperature are shown in 
Figures 2.3 through 2.8. The salient effects of the operating parameters on 
the performance of the heat pipe are discussed below. 

2. 2. 3.1 The Source Mass. 

As expected, the thermal mass of the heat source was the most important 
parameter that determined whether the heat pipe operates in a stable or unstable 
mode. All high mass runs were unconditionally stable and generally character- 
ized by a single damped temperature oscillation after a step change in power 
followed by stable operation. On the other hand, low mass runs were found to 
be characterized in general by temperature oscillations whose amplitude and 
frequency depended to various degrees on other operating parameters. 

2. 2. 3.2 The Sink Temperature. 

Under stable operating conditions, the sink temperature has a relatively 
minor effect, namely a shift in the operating temperature level of the heat 
pipe and the source. This effect can be seen in Figures 2.3 and 2.4, where 
the input heat pipe adiabatic temperature is shown as the power is increased 
in steps from 0 to 100 watts for sink conditions of -75°F and 0°F, respectively. 
The transient profiles are seen to be quite similar, except that the overall 
temperature level is 2°F to 4°F higher at 0°F sink than at -75°F sink. 

For runs where temperature oscillations were present, lower sink tempera- 
tures had a stabilizing effect. This can be seen in Figures 2.5, 2.6 and 2.7 
which show the adiabatic temperature of the input heat pipe of the VMHP operated 
at 100 watts with a small source mass, large sensor, and the condenser close- 
coupled to the sink which was held at temperatures of 0°F, -35°F, and -70°F, 
respectively. Large temperature oscillations at O^F are substantially reduced 
and eventually completely damped as the sink dropped to -35°F and -70°F, re- 
spectively. For the lowest sink condition, the onset of temperature oscilla- 
tions after a period of stable operation can be attributed to an observed drift 
upward of the sink temperature. This observation would seem to indicate that 
the oscillations were critically damped near a sink temperature of -70°F, and 
that further lowering of the sink temperature could have resulted in continuous 
stable operation at 100 watts. 
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TABLE 2-2. 

tOLQOU.i; eraml e > 

VAP0R MODULATED TRANSIENT TEST RUNS 

O 

+-» 

1 

Sensor | 
Volume ' 

Maximum 
Temp., °F 

Input Heat Pipe Adiabatic Section 

Minimum Period Mean or S.S. j Initial 

Temp., °F Minutes Temp., °F j Overshoot, °F 

Remarks 




90 "1 

80 

1 1 .4(Hal fPeriodj 

NA j 

NA 

1 NA i 

78.5 j 

78 

I Variable j 


Kg l 



fl .4(HalfPeriod| 
i NA j 

1 «A [ 

h.0(HalfPerio4 

,3 * 

1.4 " •: 


1-7 s 

,, 

[ 1.6 I 

f — 

12 . 4 (. 64 , 1 . 8 +) ) 

! 2.1 | 

{ 1-8 I 


(851 

/38l 

f /l •' 9\ 

t 

j 

/62l 

169/ 

168/ 

> 12,5/ 

i 

1681 


86.5- 79 = 7.5 No oscillations. 

88.5- 82 = 6.5 No oscillations. 

80-75 =5 No oscill . .Overshoot from 50 watts. 


Single Damped oscillation. 

No oscillations. 

Irregular oscillations, then burn out. 


No oscillations. 

....... .. , No oscillations. Overshoot from 50 via 

78-75 =3 I No oscillations. Overshoot from 75 wa 

92-84 = 8 jNo oscillations. 

84-80 =4 I No oscillations. Overshoot from 50 wa 


I 84-82 
! 81-78 


| 71-49 

| 53-35 


, 118-83 = 


111-104 


128-126 


144-138 = 6 


J No oscillations. 


Overshoot from 50 watts. 
Overshoot from 75 watts. 


Single damped oscillation. 

2 No oscillations. Voershoot from 50 watts 

3 No oscillations. Overshoot from 75 watts. 

Single damped oscillation. 

Single damped oscillation. 

Intermittent damped oscillation. 

15 Steady oscillations. 

Steady oscillations. 

14 Steady oscillations. Overshoot from 75 watts. 

15 Steady oscillations. 

Variable amplitude oscillations. 

Variable amplitude oscillations. 

23 Irregular amplitude oscillations. 

22 Steady oscillations. Overshoot from 50 watts. 

28 Intermittent stable oscill. Overshoot from 75 

21 Steady oscillations: 

! Steady oscillations. 

! /Initial damped oscillations. 1 
IResidual steady oscillations. I 

Steady oscillations. 

i 

Steady oscillations. 

Steady oscillations. 

7 No oscillations. 

No oscillations. 

No oscillations. 

2 No oscillations. 

No oscillations. 

No oscillations. 

6 No oscillations. 

No oscillations. 

No oscillations. 


NA = Not Applicable 
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2. 2. 3. 3 The Power Input. 

With a large source mass, the response of the heat pipe to a step change 
in power was characterized in general by a single damped temperature oscilla- 
tion with a period of less than 10 minutes followed by steady operation at a 
temperature level that decreased with increasing power. The decrease of the 
operating temperature at higher power is attributed to the higher operating tem- 
perature of the sensor which results in lower pressure drops across the valve. 
Small source mass runs showed that the response of the heat pipe to a step 
change in power is generally characterized by temperature oscillations. The 
amplitude and frequency of these oscillations depend on the power level and 
other operating parameters, particularly the volume of the sensor and the sink 
temperature. From Table 2.2, it can be seen that the period of the oscilla- 
tions decreases with increasing power, which tends to favor stability. The 
stabilizing effect of high power can be seen in Figure 2.8, where the input 
heat pipe adiabatic temperature is shown as the MVHP was operated with the 
larger sensor volume and the condenser close-coupled to the sink at -70°F. As 
shown, the large, variable-amplitude oscillations observed at 50 watts become 
stable at 75 watts and critically damped at 100 watts. 

Analyzing the results from Runs 19 through 24, it can be seen that high 
power can have a dampening effect on temperature oscillations provided the 
sink temperature is sufficiently low. 

2. 2. 3. 4 The Sensor Volume. 

One of the original objectives in the experimental investigation of the 
thermal response of VMHP was to determine whether or not a larger sensor volume 
could result in increased control sensitivity. The experimental results were 
only partially successful in resolving this issue. As seen in Table 2.2, the 
results with the larger sensor, Runs 25 through 42, show that the set point 
increased continuously during this phase of testing. As a result, a viable 
comparison between the two sensors cannot be made. 

Comparing Runs 13 through 24 to Runs 25 through 42, it can be seen that 
under similar operating conditions, the performance of the heat pipe with the 
smaller sensor is significantly more stable than with the larger sensor. This 
fact can be construed as resulting from the higher sensitivity of the larger 
sensor. 


14 



High source mass Runs 34 through 42 show that the set point shifts upwards 
as the sink temperature is increased. It is surmised that the larger fluid 
volume expansion in the larger sensor, particularly during source temperature 
overshoots, could have caused permanent enlargement of the volume of the bellows 
subsequently requiring higher sensor operating temperatures to open the valve. 
This seems to be supported by the experimental results. If this is correct, 
the larger sensor exceeded the design constraints of the valve/bellows assembly 
of the present vapor-modulated heat pipe. 

2.3 Investigation of the Use of Non-Condensible Gas 

During tests of the VMHP conducted under the previous contract, reduction 
of the sink temperature below the freezing point of the working fluid caused 
the heat pipe to fail. A proposed solution to the freezing-point limitation 
was to add non-condensible gas into the output heat pipe. 

2.3.1 Theoretical Analysis 

In principle, as the active portion of the output heat pipe falls in tem- 
perature, the gas expands and blocks an increasingly longer region of the con- 
denser. While the gas blocked region can freeze, the active portion is kept 
above the freezing point, and the condensing liquid is free to circulate. The 
flat-front theory can be used here to estimate the amount of gas required. 
Assuming that at maximum heat load and highest sink temperature the active 
vapor temperature of the output heat pipe is 70°F and a gas reservoir is at 60°F, 
and under minimum load and the lowest sink conditions the active vapor tempera- 
ture is -100°F with the sink at -200°F, then the required reservoi r-to-condenser 
ratio is 0.13. Because of this small ratio, the gas can be simply stored at the 
end of the condenser. For a condenser section 24 inches long, the required 
amount of gas is approximately 5 x lO - ’’ lb-mole. 

2.3.2 Experimental Results 

During the transient response test runs, it was found that some gas was 
present in the output heat pipe as indicated by temperature gradients in the 
condenser section. This gas was vented prior to the final transfer of the cal- 
culated amount of nitrogen gas into the output heat pipe. 


The VMHP was instrumented with thermocouples as in previous test runs, and 
was tested with the high source mass and the condenser thermally close-coupled 
to the cold plate which was cooled with liquid nitrogen. The performance of 
the heat pipe under various heat loads and sections of the condenser below the 
freezing temperature of the working fluid is shown in Figure 2.9. At 20 watts, 
the adiabatic temperature of the input heat pipe is seen to oscillate about a 
mean temperature of 80°F with an amplitude and frequency which remain steady 
for more than 2.5 hours with most of the condenser frozen. As the power is 
increased to 50 watts, the temperature oscillations are damped followed by 
steady operation with most of the condenser still frozen. At 100 watts, the 
gas front is seen to move past TC 13. Steady temperature control can be ob- 
served for 40 minutes with more than half of the condenser below the freezing 
point. Thereafter, the flow of liquid nitrogen was interrupted to allow the 
condenser to warm up. For the next two hours, as the sink temperature increases 
from -157°F to 25°F, the input heat pipe adiabatic temperature is seen to re- 
main essentially undisturbed by the large changes in sink conditions pointing 
to the excellent control characteristics of the VMHP. 










2.4 Theoretical Investigation of Transient Response 
2.4.1 Discussion 

Transient test results of the vapor modulated heat pipe show cyclic 
adiabatic temperature profiles that can be used as a basis for formulating 
an analytical model of the valve performance. Figure 2.10 is copy of data 
obtained during the freezing point testing. Portions of the curve are 
labeled for ease of discussion. 

The postulated phenomena causing this profile are: At point A the valve 

closes and the vapor space in the evaporator starts to be pressurized by 
the still -evaporating liquid. The rise in vapor pressure causes a rise in 
wetted-wall temperature, and the heat flow from the evaporator mass to the 
evaporator wall is reduced, while the source heat flow continues. The 
difference between the two heat flows causes the evaporator mass to be 
heated transiently. The rate of heating is reduced by the remaining evapor- 
ative cooling as the vapor is compressed and blows liquid out of the evapor- 
ator wick and vents through the wick to the condenser. 

The continued evaporation and loss of liquid blown through the valve 
bulkhead causes the wick to desaturate. With continued desaturation, 
eventually a vapor-liquid pressure difference (stress) is achieved sufficient 
to prevent capillary pumping up the evaporator wall grooves. The groove 
dry-cut at point B eliminates the bulk of the evaporative cooling, and the 
bulk of the heat flow goes into transient heating of the evaporator mass. 

The slope of the temperature-time curve becomes steep . 

As the evaporator mass heats up, the thermosat volume heats up too, but 
with a time lag. Eventually the thermostat heats up sufficiently to open 
the valve enough for the meniscus bridging it to rupture, and vapor evaporated 
on the evaporation side of the valve is vented, allowing liquid to rewet the 
wick. The rewetting of the wick gives some cooling, perhaps from boiling in 
the wick, and the cooling reduces the rate of temperature rise at point C. 

The full length of the wick is rewet at point D. 

The wick saturates sufficiently to allow the grooves to reprime, rewetting 
the evaporator wall at point E. Now the cooling is rapid. As the evaporator- 
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mass temperature falls, so does the thermostat temperature, but with a time 
lag. Eventually the thermostat closes the valve at point A', and the cycle 
repeats. 

It is desired to computer-model the transient behavior in order to 
quantify the proceeding qualitative explanation. If the computer model 
agrees reasonably well with the observed behavior, then the phenomena are 
probably understood, and the effects of various system parameters such as 
evaporator mass and thermal resistance can be determined. 

2.4.2 Mathematical Model Elements 

Equipment and Heat Source. For simplicity the equipment and heat source 

O 

will be modeled as a single lumped capacity C-| heated at rate Q and cooled by 
thermal coupling to the vapor space at temperature T v . The thermal coupling 
to T v is by a main path to the evaporator wall and an auxiliary path via the 
evaporator wick. A minor shunt conductance representing conduction to the 
condenser may be assumed as well, of course, the thermostat acts as a minor 
heat sink during transient warmup. Figure 2.11 illustrates the model and 
Eq. (1) gives the governing equation: 

dT l . . 

C 1 Ht“ = Q " Q l-e ' Q l-2 (1 ) 


where 


Q-e 


T 1 - T e 


V 


+ R + R 
o e 


*♦ 


1 


S 1 a. 1 

R_ n . 1 


S R., + 


v i_ + i_ 

R n R ,, 
g w 


( 2 ) 
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C-| Equipment-mass thermal capacity 

C£ Thermostat thermal capacity 

R. Input heat pipe resistance 

R t Thermostat coupling resistance 

Rg Grooved wall resistance (negligible) 

R,, Wick resistance 

w 

R s Shunt resistance 

R y Vapor valve resistance (nonlinear) 

R Q Output heat pipe resistance 

R g External condenser resistance 


Figure 2.11. Thermal Model 
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and 


' T 1 - V 


Thermostat . The thermostat is also modeled as a single lumped 
capacity. 


C _i _ I ( T . T ) 

c dt " R t 2 


Grooves . The grooves are modeled with a switch that is either 

on or off and resistance R . Of course, the grooves could be modeled as 

9 

a number of such subelements, but one is thought sufficient for a first- 
stage model. The switch U g is equal to either 0 or 1 and is imagined to 
be a function of vapor-to-liquid pressure difference P y - which in turn 
is a function of the liquid remaining in the wick m-j . 


°* m l < ni crit, grv ^ '"small 


1, m-j > ni cr ^^ g rv ( ' m sma i ] ) 

Valve . For simplicity the valve is modeled as either open or closed. 
When the valve is opened R v = 0. (This crude feature of the model can be 
refined later.) When the valve is closed the vapor must blow through the 
large pores in the wick, and the vapor flow is modeled as subsonic Darcy 
flow through a porous plug. 


In the t.wo-pore model to follow. 
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. T v - T b . T v - T b . 


T v- T b 


Q t. _ a 2D 


v fg h f 9 p v A ^7TT [P v ap (T v ! - P vap ( V ] 




T v - T b 


&WV - WV^ 


T„ < T 


set 


(6a) 


R v = “'V T set 


(8b) 


For convenience the Clausius - Clapeyron equation may be used for P 


vap 


(T). 


P vap' T > 


exp 





(7) 


It is necessary to find and from Figure 2. For example 

t * T i - tii-. - < T i - v/y^/'V 1 + y 1 ” (si 

Wick. The mass of liqu'i in the wick controls groove wetting through 
Eq. (5), and through a relation yet to be specified. The liquid remain- 

ing in the wick is controlled by liquid flow out (or in) and the rate of 
vaporization. A two-pore model may be used for simplicity in representing 
the wick. The wick is imagined to have a cross section of A of big pores 
of size D and one of a_ of small pores of size d. (1) The wick can oe fully 
saturated in which case both areas are filled over the entire length L. 

(2) The wick can ha big-pore de f .iterated in which case A is only partly 
filled uniformly over the entire length L. (3) The wick can be small-pore 
desaturated in which case area A_ is dry and area a is fully filled over 
only part of length L. 



The equation for 

I 


The quantity m eyap 


is simply 
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dm.| 

dt "Wap + "'in 


is given by the circuit in Figure 2.11 


"Wap Wap^fg 


T i - T b 


*evap 


R v + rtr 

R g R w 


T 1 * T 1 ' 


(9) 


(10a) 


(10b) 


(10c) 


T b " T e + $i-e (R a + R o> 

The quantity m in is more complex. Three cases corresponding to the three 
states of the wick must be considered. In general, however. 


m 


in 


'4 D 2 p 


i 4 d S 

f wet A + Ktu.L 


1 wet 



4cT - AP 
D A V 

J 

L men 


(11) 


AP = P (T ) - P (T ) 
v vap' v' r vap u b ; 


For the wick to be full, the valve must be open (but not necessarily 
vice versa). Under these conditions the following valves exist: 


wet 


= 1» L, 


wet 


= L, T v = T b so AP y = 0 


D > 
men 



m, 


full 


U = 
9 


w 


1 


( 12 ) 
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The meniscus diameter-of -curvature adjusts so that dm^/dt = 0. Recall that 
this situation was imagined to occur between E and A' in Figure 2.10. After 
the vclve closes the large pores begin to desaturate, because Tj heats up 
transiently due to the reduced value of R y , and eventually T y - T b increases 
to the extent that 


APV > 


4cr 

D 


Even before this condition, Eq. (9) will have gone to negative value of 
dm-j/dt. With m small < we have 


wet 


. "l - "'small 


m 


full 


m 


small 


L wet " D men 


= D, U = U =1 
g w 


(13) 


With continued operation at negative dm./dt, the value of m^ falls 
bel 0 Wm S«BU- 


n 'l < "'small tiet _ " b nien “ ^ 

v° vw l 

W L “ n 'f '"small 04) 

When the valve reopens AP y falls to zero, and dn^/dt becomes positive. 
Until m-j reaches m the conditions of Eq. (14) remain in force. When 

exceeds m snia -j ^ > Eq* (13) comes into play. Finally, when m-j reaches 
nifun> the conditions of Eq. (12) pertain. 

The value of R w needed is taken to be simply 

fc ' t u « (15) 

where R wq is for boiling in the full wick. 
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2.4.3 Computer Program Elements and Organization 

Main Program . The main program is visualized to read in the required 
data, initialize the variables, and then repeatedly call a subroutine to 
execute a Runge-Kutta forward marching in time. The time-temperature history 
is written as output. 

Runge . The main subroutine is thus visualized to be subroutine Runge, 
w’lich calls Subroutine Delta four times and averages the results. The main 
variables are equipment temperature Tj, thermostat temperature T 2> and mass 
of liquid in the wick m-j . 

Delta . Subroutine Delta computes the derivatives dTj/dt, dT^/dt, and 
dm^/dt using Eq. (1), Eq. (4), and Eq. (9) respectively. Functions or 
subroutines must be used to solve (given T. , T ? and M, ) for R , R , R , 

and T and T, . Subroutine Test may be used to find R„ and T versus m, , 

v b g w 1 

and Subroutine Relax may be used to find T y , T fa , and R . 

Test. Subroutine Test is visualized to find R_ and R versus m, bv 

g w 1 J 

employing Eqs. (12), (13), or (14) depending upon whether m^ = m f u n» 
m small - m l <m full’ orm l < m smal 1 ' 

Relax . Subroutine Relax must solve the nonlinear problem of finding 
Tjj given by Eq. (lOd) and T y given by Eq. (8), with Eqs. (lObc), and R y 
given by Eq. (6a) or (6b). It is thought that a relaxation procedure can 
be used or quasi 1 i neari zation . 

Pvap . Function Pvap computes the vapor pressure from Eq. (7) or a 
more accurate fit to vapor pressure data. 

Two-Pore-Model Parameter Estimations . 

Use of the two-pore model requires specifying six parameters: a. A, 
d, D, (Kx).,, and (Kt ) 2 appearing Eqs. (6a) and (11). There are only two 
main constraints (1) that the areas sum to the total void area and (2) 
that the permeability agree with the measured value. 


a + A = e v Vick 


(16) 


^IGXNal Pack to 
°F poor aoSr 


xVick = Hbj A + xtr)- 


(17) 
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Four arbitrary conditions must be imposed to complete the description. 
These are 


a = A 

(Kt)-j = (Kt ) 2 

D = 2d 
d = 26 


(18) 

(19) 

( 20 ) 
( 21 ) 


where 6 is the fiber diameter of a fibrous wick. The last condition is 
chosen in view of the fact that Eninger's theory for fibrous wicks has 


AP5 i 2 
a 2 

where we desire 

APd i o 
a ~ 0 

For consistency with the proposed model of groove dry-out occurring 
after partial wick desaturation it is necessary that j.ie equivalent groove 
pumping diameter D cr fall between D and d. 

where 


4q _ o_ R _ ^1/2 

n R * K cose 


cr 


1/2 


The quantity W^ 2 is the half groove width, and e^ 2 is the half groove angle. 
For a 5 mil fiber diameter 6, a 7 mil groove width and a 19 degree groove 
half angle we find compatibility, for 


D = 20 mils; d = 10 mils; D cr =14.8 mils. 
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2.4.4 Results and Conclusions 


A computer program was developed which incorporates the elements and 
the organization outlined in the previous section. A listing of the 
Fortran IV version of this program is given in Appendix B. Attempts to 
obtain physically realizable solutions failed due to numerical instabilities; 
as a result, the soundness of the postulated model could not be verified. 

It appears that the relaxation procedure used in the program is too simple 
in view of the strong nonlinearities in the mathematical model. A possible 
solution to the numerical problems is the linearization of the model by 
techniques such as quasilinearization. However, this effort was beyond 
the scope of the present contract. Since the two-pore model proposed is 
believed to be conceptually complete, continued program development is 
recommended. 


3. DEVELOPMENT OF VAPOR INDUCED DRYOUT DIODE HEAT PIPE 


) 

3.1 INTRODUCTION 

Heat pipes are self-contained devices that transport heat at a high 
rate between a source and a slightly cooler sink. In some applications, 
it is desirable to have the heat pipe transfer heat in only one direction; 
that is, when the source is warmer than the sink, the heat pipe would act 
in its normal high conductance mode, but when the sink temperature rises 
above the source, the heat pipe would turn off and act as an insulator. 

Thus, the heat pipe acts as a diode. An example is the cooling to cryogenic 
temperature of an infra-red detector on a satellite. A diode, heat pipe 
couples the detector to a space radiator that radiates to deep space. 

When the sun shines on the radiator, it becomes much warmer than the 
detector, and the heat pipe turns off, thus preventing the detector from 
warming up. 

Currently, there are three mechanisms that are used to achieve diode 
behavior in heat pipes: 

1. Gas Blockage - In the forward mode, a reservoir at the end of 
the condenser is used to store noncondensible gas. In the 
reverse mode, the gas blocks the region of the heat pipe at the 
source. 

Disadvantage - A very large gas reservoir is required, especially 
for cryogenic applications. 

2. Liquid Trap - In the reverse mode, liquid is condensed in an 
unwicked volume near the source and the wick dries out. 

Disadvantage - A large excess-liquid reservoir is required and 
considerable heat is transferred in the reverse 
mode before the wick dries out. 

3. Excess Liquid Blockage - In the forward mode, excess liquid is 
stored in a reservoir at the end of the condenser. In the reverse 
mode, the excess liquid fills the vapor space in the region of 
the source and prevents heat transfer. 

Disadvantage - A reservoir is required and heat is transferred in 
the reverse mode to evaporate and condense the 
excess liquid. Also, it is not well suited for a 
slab-wick heat pipe. 
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A fourth mechanism reported here for achieving diode behavior has none 
of the disadvantages of the other three mechanisms. This mechanism was 
developed under the present contract: a vapor-induced dryout diode using 

vapor modulation by means of a metal -foil reed valve. 

3.2 BASIC OPERATING PRINCIPLE 

The mechanism that was previously used successfully for passive feed- 
back temperature control, which is based on the vapor-induced dryout 
principle, is depicted in Figure 3.1. Although this design differs some- 
what from the one selected for fabrication, it has essentially the same 
basic features and is presented here to illustrate the diode mechanism 
and as a possible alternative design. 

A conventional non-arterial slab-wick heat pipe has a metal -foil reed 
valve in its transport section. In the forward mode, the reed valve is 
designed to open at a pressure differential which is small compared to the 
capillary pressure limit of the wick. Returning liquid from the evaporator 
bypasses the valve bulkhead through an annular transfer wick. In the 
reverse mode, the valve closes. When the resulting pressure that builds 
up reaches the capillary pressure limit of the wick, which occurs rapidly, 
liquid flows from the high pressure side of the valve, through the annular 
wick, to the low pressure side. The region that was the condenser in the 
forward mode dries out and the region that was the evaporator floods with 
liquid. 

3.3 DESIGN OF DIODE HEAT PIPE 

The design is shown in SK78007 in Appendix A. A 0.50-inch O.D. 
graded -porosity slab-wick heat pipe has in its transport section a 1-inch 
O.D. , 1.5-inch long cylinder which supports the valve bulkhead. The 
density radiation of the slab-wick with respect to length is shown in 
Figure A-l in Appendix A. 

In addition to its high transport capacity, the graded porosity wick 
has certain features that are ideally suited for the diode mechanism. The 
pressure available to open the valve is the high capillary pressure 
generated by the relatively low porosity of the evaporator region, while 
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the pressure required to induce dryout is the low capillary pressure of 
the relatively high porosity of the condenser region. 

All components were made out 304 stainless steel. The criterion for 
selecting a single material to build the heat pipe diode was the ready 
availability of the materials and manufacturing cost. The use of aluminum 
tubes connected to the valve section by means of transition tubing would 
have been desirable from weight consideration, however, the higher manu- 
facturing cost could not be justified at this stage when the main objective 
was to test the practicality of the diode mechanism. 

3.3.1 The Valve Bulkhead Design 

The valve bulkhead assembly before it was assembled is shown in 
Figure 3.2. The bulkhead has a central solid core with five orifices, 
a wick feed-thru orifice and four smaller orifices for the vapor flow 
tubes which are pressure fitted into holes in the bulkhead and pinned in 
place. The ends of two vapor flow tubes provide the seat for each of the 
two reed valves. In order to fit the valve assembly into a 1-inch O.D. 
cylinder and still be able to use relatively long foil flaps, the ends of 
the vapor flow tubes were cut at an angle placing the seat of each valve 
on a plane approximately 30° with respect to the axis of the heat pipe. 

On this plane, the vapor flow crossectional area at the valve seat is 
about twice that at the inflow, thus only small deflections of the foil 
reeds are required to fully open the valve. The vapor flow tubes have a 
thin wall which is necessary in order for the valve seat to have a small 
contact area, otherwise surface tension forces due to liquid between the 
seat and the foil would generate an excessive force tending to keep the 
valve shut. 

3.3.2 The Wick Feed-Thru 

The wick feed-thru is shown in SK78007 and also in a photograph in 
Figure 3.3. Its configuration is similar to that used in the vapor- 
modulated heat pipe. The wick on either side of the bulkhead is separated 
at the bulkhead by two layers of 250-inesh screen, which form the capillary 
barrier. The screen is held in the center of the wide feed-thru orifice 
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by two sleeves that were inserted halfway into the orifice from both sides 
of the bulkhead. The sleeves were pressed together to clamp the screen 
in place and then were tack welded to the bulkhead. The ends of the wicks 
were inserted into sleeves, trimmed flush, inserted into the bulkhead 
sleeves, pressed against the capillary barrier, and pinned in place. 

3.3.3 Reed Valve Assembly Tests 

The valve assembly was initially tested in a glass tube which allowed 
visual observation of its behavior. A Teflon bushing between the glass 
tube wall and the bulkhead provided an air-tight seal. The glass tube 
was instrumented to allow pressure and gas flow measurements. 

Reed valves manufactured from 0.0005-inch and 0.00025-inch 304 
stainless steel foils were tested, and the latter chosen due to their 
greater sensitivity to pressure differentials. 

The effective pore size of the capillary barrier in the test assembly 
was determined by a bubble-point test. The wick was saturated with acetone 
and the reed side of the test chamber was pressurized simulating diode 
behavior. The pressure continued to rise until the capillary barrier 
failed and bubbles were observed emerging from the wick. The pressure 
at failure corresponds to an effective capillary pore size of 0.0035-inches. 
This pore size is sufficiently small to ensure that the abnormal evaporator 
will dryout when the valve closes in the diode mode. 

Preliminary tests were performed in the forward mode in which the 
wick was saturated with acetone, the valve seats were wet, and helium 
gas was used to simulate the vapor flow. It was observed that only one 
of the valves would open. This behavior was attributed to the permanent 
set of one of the foil reeds opposing outward deflection of the foil and/or 
large capillary forces keeping the foil adhered to the valve seat. The 
foil flaps were subsequently placed between two metal blocks and annealed. 
This operation resulted in flat foil material. 

The valve seat contact area was reduced by chemically etching the 
front section of the vapor flow tubes. 
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The valve assembly was then reassembled and subjected to further 
tests. Under conditions similar to those described above, it was observed 
that the pressure continued to rise until one of the valves opened. The 
pressure at which the valve opened corresponds to an effective capillary 
pore size of 0.026-inch. This pore size is sufficiently large to ensure 
that the liquid in the normal evaporator will not flow into the condenser 
side before the valve opens. 

For a wide range of gas flow rates, the foil was found to deflect 
smoothly without flapping. By diverting the gas flow into the other side 
of the glass chamber, the valve was tested in the diode mode. The 
pressure would rise until the capillary limit of the barrier was reached. 
Simulation of heat pipe operation in the forward and diode modes in this 
way was repeated at least a dozen times with reproducible results. 

3.3.4 Heat Pipe Diode Test Results 

The heat pipe was assembled following the procedure outlined in 
SK78007, page 4 of 4, a photograph is given in Figure 3.4. 

After a 16-hr refluxing operation at 80°C, the heat pipe was vacuum 
baked and charged with 32.6 grams of ammonia. Subsequently the heat pipe 
was instrumented and tested in the forward mode. At a 0.2-inch tilt, 
it held 135 watts and burned out at 150 watts. This was well below the 
predicted capacity of over 225 watts. In addition a temperature drop 
was observed across the va^e, which increased with increasing heat input. 
These two facts were evidence of an impaired function of either the wick 
or the valve or both. 

Since structural damage to the valve was suspected, prior to con- 
ducting a detailed tilt capacity test in the forward mode, the heat pipe 
was tested in the reverse mode. Repeated attempts to close the valve by 
applying increasing power failed inasmuch as the evaporator showed no signs 
of dryout even at 100 watts. In fact the heat pipe operated much as in 
the forward mode. 

The results clearly indicated the valve was damaged. 
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Figure 3.4 Diode Heat Pipe After Assembly 

The heat pipe was subsequently dismantled. A photograph of the valve 
assembly is shown in Figure 3.5 where it can be readily seen that the 
reed valves had been severely damaged. 

The foils were wrinkled and showed evidence of having been forced 
into the tube sections. This damage would not allow the foil to seat and 
therefore prevented the foils from stopping vapor flow in the reverse 
direction. 

One additional item of damage sustainea during manufacture was a 
hole in the wick about one inch from the bulkhead on the condenser side 
of the heat pipe. This hole was approximately 1/8" diameter extending 
approximately 2/3 of the way through the wick. The hole was darkened and 
was obviously the result of a burn. Apparently a molten drop of metal 
tell into the wick at some point during welding-most likely during one 
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of the repair welds made to close the leak in the valve section. The 
burn spot corresponded to the location of the leak. This wick damage is 
the likely cause of the relatively low capacity of the heat pipe (135 watts 
actual vs 225 watts predicted). 

3.4 CONCLUSIONS 

The heat pipe diode components were tested under simulated normal 
operating conditions and found to be a sound working mechanism. The 
assembled heat pipe diode failed due structural damage incurred under 
conditions not yet fully understood. 

It is quite likely that the damage occurred during some point in the 
processing stages. During the routine leak check operation, a leak had 
been detected in the area where the 1/2 inch O.D. tube is welded to the 
end cap of the valve section. In the process of finding the location of 
the pinhole, there was an accidental pulse pressurization of the heat pipe. 

L 
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The sudden pressure differential sustained across the valve would 
have caused the valve to bend open, but a resultant shock wave could have 
impelled the foil back against the valve seat with a force sufficient to 
cause the observed damage. 

Another possibility is that damage occurred during the refluxing 
step during which the pressure fluctuations associated with the unstable 
nature of pool boiling could have caused the foil to implode against the 
valve seat. 

The successful testing at the component level of the metal -foil 
reed valve and the high capacity of the metal-fiber wick system still 
make this diode mechanism attractive. Re-work of the parts and careful 
processing of the diode seem warranted to more fully evaluate its 
performance. 


} 
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APPENDIX A 

CONTENT 

SK78007 - GRADED POROSITY HEAT PIPE DIODE ASSEMBLY 

SK741001 - TUBE, GROOVED 

SK741024 - END CAP AND FILL TUBE 

SK75061 - VAPOR-MODULATED HEAT PIPE INSTRUMENTATION 

FIGURE A-l - HEAT PIPE DIODE GRADED-POROSITY WICK 
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VOLUME DENSITY OF WICK 




APPENDIX B 

Listing of computer program VMHP to 
simulate vapor-modulated Heat Pipe 
Performance 


0010Q 

00110 

00120 

00130 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

C0370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00490 

00500 

00510 

00520 

00530 

G0540 

00550 

00560 

00570 

00580 

00590 

00600 

00610 

00620 

00630 

00640 


PR06RAH VMHP< INPUT, OUTPUT, TAPE5=INPUT,TAPE6=0UTPUT) 
C0MM0N/ALL/C1,C2,RI,RTS,RG,RU,RS,RV,R0,RE f T1,T2,TE,TB,TV,TI, 
11 , LUET'LV^HEM'DMfAA.HL'HLSHAL'HLFULMG'UMIE, 812,0, 
2VISV, VISL,RHOL,RHOV,HFG, SIGMA, KTAO,TIME,DT 
3,»PV,TSET,UML 

REAL L ,LUET ,LV, ML ,MLSMAL, MLFULL ,KTAO 
M«60 

CALL BATA 
TIME=~DT 

100 00 101 1-1 , N 
TIME=TIHE+DT 
CALL RUNGEdl ,T2,ML,DT) 

CALL OUTPUT 

101 CONTINUE 
STOP 
END 

C 

C UPDATE T1,T2,HL 
C 

SUBROUTINE RUMGE ( Y10,Y20,Y30,H> 

CALL DELTA < Y10,Y20,Y30,DYf ,DY2,DY3) 

SUM1 =DY1 
SUM2=BY2 
SUH3=DY3 
Y 1 =Y1 0+DY1 *H/2. 

Y2=Y20+DY2*H/2. 

Y3-Y30+DY3*H/2. 

CALL »ELTA(Y1 ,Y2,Y3,DY1 ,BY2,DY3) 

SUM I =SUH1+2.-* : DY 1 
SUM2=SUM2+2.*DY2 
SUM3~3UH3+2.*DY3 
Y1=Y10+BY1*H/2. 

Y2-Y20+DY2+H/2. 

Y3=Y3D+DY3*H/2, 

CALL DELTA* Y1 ,Y2,Y3,DY1 ,DY2,DY3) 

SUH1=SUM1 +2. *»Y1 
SUH2=SUM2+2.*DY2 
00480 Y 1 =Y1 0+DY1 *H 

Y2=Y20+BY2*H 
Y3=Y33+DY3+H 

CALL DELTA ( Y1 ,Y2,Y3,DY1 ,BY2.DY3> 0 

Y10-YI0+(SUM1+DY1 )*H/6. 0^ V 

Y20=Y20+< SUM2+DY2 )*H/6. 

Y3Q=Y30+<SUi13+BY3)4H/6. 

RETURN 
END 
C 

C CALCULATE DTi/HT,BT2/l!T, AND BML/DT 
C 

SUBROUTINE BELT A( XI ,X2,X3,DT1DT,BT2BT,DMLDT) 

COMMON/ ALL/C 1 ,C2,RI ,RTS,RG,RU,RS,RV,R0,RE,T1 , T2,TE,TB,TV,TI , 
lL,LUEr,LV,UHEN,D,BB,A f AA,NL,MLSNAL,MLFULL,UG,UU f Q1E,Q12 f G, 
2VISV, VISL,RHOL,RHOV,HFG, SIGMA, KTAO, TIME, BT 
3,DPV,TSET,UML 
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00650 

00660 

00670 

00680 

00690 

00700 

00710 

00720 

00730 

00740 

00750 

00760 

00770 

00780 

00790 

00800 

00810 

00820 100 

00830 

00340 

00850 C 

00860 C 

C0870 C 

00880 

00890 

00900 

00910 

00920 

00930 

00940' 

00950 

00960 

00970 

00980 

00990 

01000 

01010 100 
01020 
01030 
01040 
C1050 
01060 
-U070 
01080 50 

01090 
01100 
omo 
01120 
01130 
01140 
01150 
01160 
01170 


REAL L,LUET,LV,HL,HLSMAL,NLFULL,KTAO 

T1=X1 

T2=X2 

HL=X3 

CALL RELAX 

DTIDT=(Q-Q1E-Q12)/C1 
BT2DT=Q12/C2 
IFiUKL.EQ.O. ) GO TO 100 
X=UG/RG+UU/RU 
X=RO+l./X 

ENEVAP=<TI-TB)/X/HFG 

C01=4.*RHQL/KTA0/VISL 

CQ1=CQ1*(IiB**2*FUET*AA/L+Ii**2*A/LUET> 

C02=4.*SIGHA/HMEN 

EHIN=C01 *(C02-BPO)*4. 1 7E8 

I'HLDT=-EHEVAP+EHIN 

RETURN 

DHLDT=0. 

RETURN 

END 


CALCULATE CUE, 01 2,0, AND RV 


SUBROUTINE RELAX 


COHMON/ALL/C1,C2,RI,RTS,RMU,RS,RU,RO,RE,T1 , T2, TE, TB , TV, TI, 
1 L, LUETjLV, DMENjD, DD, A, AA, HL , HLSHAL,HLFULL,UG,UU,Q1 E,Q1 2,0, 
2VISV, VISL, RHOL,RHOV,HFG,SIGNA f KTA0,TIHE r DT 
3,DPV,TSET,UHL 

REAL L,LUET ,Ly,ML ,HLSMAL,HLFULL ,KTA0 

CALL TEST 

Q12=(T1-T2)/RTS 

IFCT2.GE.TSET) GO TO 200 

C01=2 U *RHO»*BFG:«-AA*PI'**2./KTAO/LV/VISV 


TV=T0+5. 

TB--TB-5. 
icouh r=o 

TyOLU=TU 
ICQUNT=IC0UNT+1 
IFCICOUNT.GT. 10) STOP 
TBQLB-TB 
PVTV-PV(TO) 

PVTB=py(Tn) 

iipy-pyiy-poTB 

RV- ( T V-TB ) / ( C0 1 *DP 0*4 . 1 7E8J 

RT=UG/RG+UU/RU 

RT=1 ./RT+RV 

RT=1 ./RS+1 ./RT 

RT“1 . /RT+RI+RO+RE 








Q1E-(T1-TE)/RT 
TB=T£hT1 £* CRE+RO) 

TI=T1~Q1E*RI 

TV=TI- ( 01 E-C TI-TB) / Kti) /( UU/RU+IIG/RG) 
ERTB=ABS ( 1 . - C TBHfcO. )/( TB0LD+460 . ) ) 
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01180 

01190 

01200 

01210 

01220 

01230 200 

01240 

01250 

01260 

01270 

01280 

01290 

01300 101 

01310 

01320 C 

01330 C 

01340 C 

01350 

01360 

01370 

01380 

01390 

01400 

OHIO 

01420 

01430 

01440 

01450 

01460 

01470 • 

01480 

01490 

01500 101 

01510 

01520 

01530 

01540 

01550 

01560 100 

01570 
01580 
01590 
01600 
01610 
01620 
01630 
01640 C 
01650 C 
01660 C 
01670 
01630 
01690 
01700 


ERTV=ABS ( 1 . - < TV+460 . ) / ( T VQLD+460. ) ) 

IFIERTB.LT.. 001. AND.ERTV.LT. .001) GO TO 101 
TB=(TB+TB0LB)/2. 

TV=(TV+TV0LD)/2. 

GO TO 100 

RT=UG/RG+UU/RU 

RT=RT+1 ./RS 

RT=1 ./RT*RI*RQ+RE 

G1E=(T1-TE)/RT 

TB=TE+Q1E*<RE+R0) 

TI=T1-Q1E*RI 

TV-TI~(G1 E-(TI-TB)/RS)/(UU/RU+UG/RG) 

RETURN 

END 

CALCULATE UG,,UU f UHL,DNEN,FUET, AND LUET 
SUBROUTINE TEST 

COHMQN/ALL/CI ,C2,RI,RTS,RG,RU,RS,RV,R0,RE,T1 ,T2 f TE, TB,TV,TI , 
1L,LUET f LV,DHEN,D f HD,A,AA,NL,HLSNAL,HLFULL,UG,UU t Q1E,QI2 f Q f 
2VISV,VISL,RH0L,RH0y,HFG,SIGHA,KTA0,TINE,HT 
3,HPV,TSET,UNL 

REAL L, LUET, L V, ML, HLSMAL , NLFIILL, KTAO 
UNL=1 . 

IF(HL.EQ.MLFULL) GO TO 100 
IF(riL.GT.HLSHAL) GO TO 101 
FUET=0. 

BMEN-D 
UG~0 . 

UU=HL/MLSHAL 

LUET~UU*L 

RETURN 

FUET-(HL-HLSHAL)/<MLFULL--HLSMAL) 

LUET=L 
BHEN-IiB 
UG-1 , 

UU=v1 . 

RETURN 

IFtT2.GE.TSET) UHL-0 - 
ONE N= Mi 
FUET-1 . 

LLET-L 
UG-1 A 
UU=1 . 

RETURN 

END 

CALCULATE VAPOR PRESSURE 


FUNCTION PViZ) 

BATA A I ,A2,A3,A4,A5 

1/13. ? 2374 , -4. 92 17-4 , _ . 0650 I SE- 1 , -7 . 579597E-2, 0. / 
X=Z+4<‘.0. 
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01710 
01720 
01730 
01740 
01750 
01760 
01770 
01780 
01790 
01800 C 
01810 C 
01820 C 
01830 
01840 
01850 
01860 
01870 
01880 
01890 
01900 1 

01910 
01920 
01930 2 

01940 
01950 
01960 C 
01970 C 
01980 C 
01990 
02000 • 
02010 
02020 
02030 
02040 
02050 
02060 
02070 
02080 
02090 
02100 
02110 
02120 
02130 
02140 
02150 
C2160 
02170 
02180 


X=1000./X 
X2=X*X 
X3=X2*X 
X4=X3*X ' 

Y=A1+A2*X+A3*X2*A4*X3*A5*X4 
Y=EXP( Y) 

PV— 1 44 - * Y 

RETURN 

END 

PRINT OUTPUT 




SUBROUTINE OUTPUT 

COMMON/ ALL/C1 ,C2,RI,RTS,RG,f^i RS,PV,RD f SE,T1 ,T2,TE, TB,TV r TI , 
1 L , LU E T , L V , D HE NyD , D D , A , A A , r'L v Hi S M ki, * ULL , UG,UU,Q1E, 012,0, 
2VISU, VISL,RHQL,RHOU ,HFG,SiGHA, KTAO,tIME,DT 
3,DPV,TSET, UHL 

REAL L,LUET,LV,NL,HLSHAL,MLFULL,KTAO 
IFtTIME.ED.O.) PRINT 1 

FORMAT ( 1 X, *TIME (SEC )*,6X f *TBLOCK(F)^ « t *'fS£feSO»?F I*. 
15X,*TVAP0R(F)*) 

PRINT 2, TIME, T1 , T2 ,TM 
FORMAT ' 1 X ,4 (El ! .4 , 4X) ) 

RETURN 

END 


INPUT DATA 
SUBROUTINE DATA 

CDMH0N/ALL/C1 ,C2,IU,RTS,RG,RU,RS,RV,R0,RE, T1 ,T2,TE,TB, TV,TI, 

1 L,LUEf ,LV,DMEN,D,3D,A, AA,HL,HLSHAL,HLFULL,UG f UU,G1E,Q12,Q, 
2UISV,VISL, RHOL , RHQV, HFG, SIGHA, KTAO ,TIHE, DT 
3,DF'V,TSET,UML 

REAL L,LUET f LV,HL,MLSHAL,HLFULL,KTAO 
ACCEPT DT , TE 

DATA 51 , C2, h 1 S, RI , K5, RU , RG, RD , RE 

1/,5?,6.00E -3,2. 1 0, 4. 40E-?, 115. ,3.E-1 ,2.20E-2,2.90E-2,2.90E-2/ 

DATA A, AA,H, DD, L , LV , MLSMAL , MLFULL , KTAO j 

l/5.35£-4,5.35E-4,6.95E-5,1 .39E-4, .5,7.5E-2, 1 .07E-2,2. 1 4E-2, 1 52./ ' 
DATA RHOL, RHOV, SIGflA, V ISO, VISL,HFG 
1 /37. , , 705, 1 . 1 8E-3,2 , 1 5E-2,2.83E-I ,478./ 

DATA T 1 , T2 , T SET , Q 
1/75, ,75. ,75. ,341,/ 

HL-nLFULL 

TU=T1 

TB--TE 

RETURN 

END 
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